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Abstract. The g-sum x © q y = x + y + (1 — q) xy (x (Bi y = x + y) and the q- 
product x ® q y = [x 1-9 + y x ~ q — 1]~ (x ®\ y — xy) emerge naturally within 
nonextensive statistical mechanics. We show here how they lead to two-parameter 
(namely, q and q') generalizations of the logarithmic and exponential functions (noted 
respectively hx q>q i x and e^ ql ), as well as of the Boltzmann-Gibbs-Shannon entropy 
Sbgs = — k Yli=i Pi ^ n Pi (noted Sg <g '). The remarkable properties of the (q,q')- 
generalized logarithmic function make the entropic form S q ^ q i = kY2i=iPi^ n q,q'0-/Pi) 
to satisfy, for large regions of (q,q'), important properties such as expansibility, 
concavity and Lesche-stability, but not necessarily composability. 



PACS numbers: 02.10.-v, 02.70.Rr, 05.90.+m 
1. INTRODUCTION 

Boltzmann-Gibbs (BG) statistical mechanics provides a powerful tool for understanding 
how fast microscopic physics with short-range interactions affects the thermodynamics 
of a system on much larger space-time scales. It is based on the (Boltzmann-Gibbs or 
Boltzmann-Gibbs-Shannon) entropy 

W W - w 

Sbgs = -k^Pi^Pi = k^Pi^— (^^ = 1 )' ( 1 ) 

t=l i=l ^ l i=l 

hence, if Pi = l/W(Vi), S BG s = k\nW. 

However, the BG theory turns out to be inadequate for various complex natural, 
artificial and social systems. Such is, for instance, the case where a zero maximal 
Lyapunov exponent is present (e.g., long-range-interaction many-body Hamiltonian 
systems, nonlinear dynamical systems at the edge of chaos, and others). Typically, 
such situations are dominated by power- laws instead of exponential distributions [H12]- 
In order to deal with at least part of such systems, it was proposed in 1988 [3J, a 
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generalization of the BG theory, referred to as nonextensive statistical mechanics [HE]. 
It is based on the entropy 

1 _ q w 1 w 

Sq = k — ^fPi = kJ2p> q (-) (gerc,5> = i), (2) 
q 1 i=\ Pi i=i 

(hence, if Pi = 1/W (Vi), S q = k ln g P^) where 

x i-<? _ i 

lia q x = — (hix X = In re) . (3) 

This entropy succeeds to widen the range of applicability of statistical mechanical 
concepts to many systems. It fulfills thermodynamically relevant properties such as 
expansibility, composability, Lesche-stability and concavity (for q > 0). The roles played 
in the traditional theory by the logarithmic and exponential functions are played, in the 
generalized theory, by the so-called q- logarithm, as given by Eq. (3), and its inverse, 
the g-exponential [6]: 

ej = [1 + (1 - q)x)f« (ef = e«), ( 4 ) 

where [•••]+ vanishes if its argument is nonpositive. This naturally led to the construction 
of a family of functions based on the g-generalized logarithm and its inverse function [7j. 

Consistently, nonlinear generalized algebraic forms emerge (see also [8]), namely 
the q-sum [9] and the q-product [9j [10] respectively defined as follows: 

x @ q y = x + y + (1 — q)xy (x® 1 y = x + y), (5) 

and 

x(g) q y= (x 1 ^ + y 1 ^ - l) ~ (x (8>i y = xy) . (6) 

Additionally, with these mathematical generalizations of the sum and the product, a 
new algebra could be constructed which appears to have some peculiar properties not 
found for other algebras [HI HD]- The g-product and the g-sum maintain the structure 
given by Eq. (JTj) by converting the sum and the product into these non-linear operators. 
We can straightforwardly verify the following properties: 

ln(x • y) = \nx + Iny (7) 
hiq(x -y) = ln q x® q \n g y (8) 
ln q (x ® q y) = \n q x + \n q y (9) 

In the following section (Sec. 2) we show that an uniquely further generalized logarithmic 
function, denoted ln g ^ x, exists which enables the unification of Eqs. ([319]) as follows: 

ln q>q/ (x ®g y) = ln 9i9 / x ® q > ln 9i9 / y . (10) 

In Sec. 3 we present some general properties of this function. In Sec. H]we introduce a 
two-parameter generalization of the BGS entropy, and discuss various of its properties. 
We finally conclude in Sec. 5. 
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2. DERIVATION OF THE TWO-PARAMETER GENERALIZATION OF 
THE LOGARITHMIC FUNCTION 

Eq. ([TO]) constitutes a strongly nonlinear functional equation. A general analytic 
expression of lm.^/ x needs to satisfy the limiting cases ln^i x = hii i(? x = ln ? x. 

In order to simplify Eq. ffTUj) . we set x = y and, inspired by the expression on the 
left hand side of Eq. flU]), write g(z) = g(\n q x) = ln 9i9 / x. By inserting g(z) in Eq. ({TO]) , 
this functional equation transforms into the much simpler form 

g(2z) = 2g(z) + (l-q')[g(z)} 2 . (11) 

We then verify that the ansatz g(z) = j^(b z — 1) solves Eq. ffTTl) exactly, thus leading 
to, 

ln q , q ,x = ^ (b^i^- 1 ) - l) , (12) 

where b corresponds to a constant that still needs to be determined. It is important to 
point out that \n qq '(x) solves Eq. (ITU]) also for the more general case x ^ y. 

In order to determine b, we focus on the point ln q A i 1 = 0, which is expected 
to hold for all b and all (q,q f ), analogously to the fact that ln 9 1 = for all q. We 
assume also that, at this point, the entire family ln g q i x has one and the same slope, 
i.e., -jj ln ?)9 ' x\ x =i = 1, once again in full analogy with ln 9 x. Then, we verify that our 
new function, Eq. (fT2|) . takes the following more specific form: 

x = 

Indeed, we can prove that ln 9i9 / x provides a generalized function of the g-logarithm with 
b = exp(l — q') by setting either q or q' equal to unity. Eq. f|T3|) can be rewritten as 
follows: 

ln q ^x = ln q/ e ln " x . (14) 

The relations 

ln g l x = lnx j(? x = hig x (15) 

are easily recovered by evaluating Eq. f|T3|) in the limits q — > 1 and q' — > 1. 

In Fig.fUwe illustrate the two-parameter generalized logarithmic function for typical 
q = q' cases. Typical q ^ q' cases are illustrated in Fig. [2l 

3. PROPERTIES 

In this section we present some general properties of the two-parameter generalized 
logarithm, some of which will be useful in the discussion of Sec. HI 

As can be obtained from Eq. (TT3"]) (and observed in Figs. [1] and [2]), hi qq ix 
monotonically increases with x. Thus, an unique inverse function of ln g <q i x, namely 



exp 



1-q 1 



(x 1 ^ - 1) 



(13) 
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the two-parameter generalized exponential function, exists for arbitrary values of q and 
q>: 

+ ^ Ml + (1 - q')y}] ^ = er q ' V • (16) 

There are two cut-offs, namely for y(l — q') > 1 and for {...} < 0. This can be easily 
understood looking at the shape of the generalized logarithm in Figs. [1] and [2j ln M ' x 
saturates for x — > oo if at least one of the parameters is greater than unity (either q > 1 
or q' > 1). Moreover, it also approaches a finite value for x — > if at least one of the 
parameters is less than unity (either q < 1 or q' < 1). 

It is remarkable that the functional form of the present two-parameter 
generalization of the logarithmic function contains the standard exponential function, 
e x . The reciprocal happens with the functional form of its inverse function, exp 9 g /(x), 
which contains hix. 

Many functions, e.g. sinx and cosx, can be expressed in terms of the exponential 
function. Similarly, it should be possible to construct a new family of functions using 
Eqs. (p|) and (fl6]). 

Let us calculate here the first derivative of ln q A i x: 



exp -— *■ (x l -i - 1) . (17) 



d\n. qA ,x fl-q' 

- = x q exp 

ax \ 1 — q 

The slope of ln q>q / x remains positive for all x > 0, which is consistent with the already 

mentioned monotonicity of the generalized logarithm for all values of (q, q'). 

The shape of hi^g/ x in Figs. [1] and [2] suggests that there may exist some symmetry 

with respect to the point (x,y) = (1,0). Indeed, the transformation q — > 2 — q, well 

known to be an important relation in nonextensive thermostatistics [1], makes ln q q rx 

to satisfy the following property: 

ln 9i9 / - = - ln 2 - qt2 -q' x. (18) 

4. TWO-PARAMETER GENERALIZATION OF THE 
BOLTZMANN-GIBBS-SHANNON ENTROPY 

We construct here an entropic functional based on the two-parameter generalization of 
the standard logarithm. 

Along the lines of Eqs. (1) and (2) we can construct a new entropic functional as 
follows: 

w w 

(19) 



w n , vv / , 

S q , q/ = £;5>hw - = T^-,H Pi (e^^ 1 ^ ~ 1 ) 



(with Y^LiPi = 1) hence, if pi = 1/W (Vi), S q ^ = k\n qs , W. 



Lesche- stability (or experimental robustness): The functional form of \n q g / x being 
analytic in x, this entropy is Lesche-stable (experimentally robust) [11] . 
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Figure 1. Illustration of the generalized logarithm, Eq. (Q2|), setting q = q' , in linear 
scales (left) and semi-logarithmic scales (right). 



Expansibility: A zero-probability (pi = 0) state does not contribute to the entropy if 
at least one of the parameters is greater than unity (q > 1 or q' > 1). Note that q > 
(q 1 > 0) in the limit q' = 1 (q = 1). Therefore S m >(p 1 ,p 2 , ...,p w ,0) = S q>q >(p 1 ,p 2 , --,Pw), 
i.e., S qtq i is expansible within the above given limits for the two parameters. 

Concavity: An important requirement for an entropic function is that its curvature 
does not change with {p{\ (sufficient although not necessary condition for being concave 
or convex, i.e., an entropy can still have a fixed sign of its curvature for parameter values 
not satisfying this condition). Figs. [3JA.-D illustrate S q>q i for a system of two states with 
probabilities p and (1— p) for typical values of (g, q'). In Fig.[3A (Fig. [3]B), the curvature 
of S q>q ' is negative (positive) for all (q, q'), and the extremum lies at p = 0.5 as expected. 

In order to determine the regions of the parameter space where the curvature of 
S q>q i does not change, we focus on the curvature of Pi \a. q>q i{l/pi). We check whether 
there is any change of the sign of the curvature in the interval < Pi < 1. We verify 

|U 1 = e^C*- 1 " 1 ) [Pl q ~\l ~ q') - qpl 2 ] ■ (20) 
ap i Pi 

Concavity is guaranteed if Eq. ( |20l) remains negative. This is fulfilled if q + q' > 1 in 
the limit Pi — 1, and q' > 1 if q < 1 in the limit Pi — > 0. Analogously, convexity is 
guaranteed if q + q' < 1 in the limit pi = 1, and q,q' < 1 in the limit p« — > 0. 

Fig. H] exhibits the various regions for the curvature with respect to q and q' . We 
can see that two regions (III and IV) may have non fixed curvature sign. An analysis of 
the shape of S q ^ q i with W = 2 for several values of q and q' confirms this but suggests 
that parts of regions III and IV may still have fixed curvature sign (see also Fig. [3] for 
examples of S qtq > corresponding to these regions). 

About compos ability: Consider a system A + B composed by two probabilistically 
independent subsystems A and B (i.e., pf^ B = pfpf)- An entropic functional S is said 
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Figure 2. Illustration of the generalized logarithm for typical values of the case q^ q' '. 
We set q = 0.3 and vary q' (upper left), then set q' = 0.3 and vary q (upper right). 
The lower figures correspond to fixed q = 1.7 (left) and q' — 1.7 (right). Note, that 
the curves in the left figures and the ones in the right figures can be transformed into 
each other by applying Eq. (Tl8|) . 



to be composable if its value S(A + B) depends from the probabilities {pf } {pf} only 
through S(A), S(B) and generic indices, i.e., if S(A + B) = &[S(A), S(B), indices], 
where $ is some univalued function. Any additive entropic functional (such as the BGS 
and the Renyi ones) trivially is composable (e.g., Sbgs(A + B) = Sbgs(A) + Sbgs(B)). 
The nonadditive entropy S q (q ^ 1) also is composable since ^ — = -f - © g -f- = Sg( jf^ + 
^ + (1 - q)^^- Composability is an important property for thermodynamical 
applications, not so for applications such as signal or image processing, cybernetics, and 
others. 

Let us now check whether S q ^ is composable. We start by considering a micro- 
canonical ensemble, i.e., equal probabilities. We have S^, B = k\n qtq >(WAWB), S£ q , = 
kln q>q i Wa and Si 3 , = kln qq i Wb- It straightforwardly follows that 



Two-parameter generalization of the logarithm and entropy 



7 





1.5 



t>3 



W 

5? 



0.5 



! , 


q=0.2, q'=0.9 


: \ 


q=0.4, q=0.9 


1 \ 
1 \ \ 


— - 9=0.9, ?'=0.9-i 


: ; \ ^ 
3 \ \ 


-- - ?=0.6, ?'=0.4 ! 


^ \ N v 


- - q-0.9, q>-0.4{: 






1 





0.2 



0.4 



0.6 



0.S 



Figure 3. Illustration of the two-state generalized entropy S q ^ q > for typical values 
of (q,q'). A: concavity (region I in Fig. |4|). B: convexity (region II in Fig. [4j. C: 
curvature can change (region III in Fig. 2J. D: curvature can change (region IV in 
Fig. HI). 



1 + (1 - 4 



qB 



- 1 



(21) 



The meaning of such a relation (pseudo-additivity) is uneasy to interpret. Nevertheless, 
we can present Eq. (j2ip in a simpler manner by introducing the transformation 
Y(S) = \nh + ^\n[l + (1 - q')(S/k)]\. We then obtain 

(22) 



Y(S A+B ) =Y(S A )+Y(S B ). 

Composability of the present entropic form S qyq t requires that the same relation in 
Eq. (1220 remains valid for other thermostatistical statistics such as, for instance, that 
of the canonical ensemble. However, our studies indicate that this is not the case for 
arbitrary distributions {pi}- Numerical verifications confirmed this result. Therefore 
S q>q i is in general not composable. 
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Figure 4. The different regions of the entropy with respect to its curvature 

(see also Fig. [3]). The point q = q' = 1 corresponds to the Boltzmann-Gibbs-Shannon 
entropy. The lines q — 1 and q' — 1 correspond to the nonadditive entropy 5 g . This 
diagram exhibits {q,q') — > (9', 9) symmetry although generically ^ 5V,ij- 

5. CONCLUSION 

We found an explicit generalized logarithmic function which nicely closes the picture 
obtained with the g-logarithm using the g-sum and the g-product as basic mathematical 
operations. It is possible to use the new logarithm, ln 9i? / x, and its inverse function, 
e£ in order to construct a new family of functions such as sm q>q 'X and cos m iX. It 
remains to be tested if these new findings can be useful to describe, in some way, 
some behavior of natural, artificial or social complex systems, or if they can be useful in 
applied areas such as signal and image processing, cybernetics, optimization algorithms, 
information theory, among others. Even if they do not immediately exhibit applications, 
they remain as elegant mathematical constructions. Does ln 9)(? > x constitute a "super- 
generalized" logarithm possibly useful as the basis for a more general entropy? We 
engaged this question by showing that various important properties (concavity, Lesche- 
stability, expansibility) presumably needed for a physical entropy are fulfilled in this 
case providing that some minimal restrictions are imposed to (q,q'). However, one 
crucial property failed to be valid in general, namely composability. The fact that this 
entropy is not composable leads to the assumption that it may be not able to play a 
thermodynamical role comparable to the nonadditive entropy S q on which nonextensive 
statistical mechanics is based. 
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As a last remark, let us mention that, along the present lines, Renyi' s entropy 
gR = n z^i=iP l can a j go ^ e (trivially) generalized into a two-parameter one, namely 

1-q 
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